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Abstract 

We report a detailed modeling of a mass profile of a 2 = 0.2999 massive elliptical galaxy using 30 
milli-arcsecond resolution 1-mm Atacama Large Millimeter/submillimeter Array (ALMA) images 
of the galaxy-galaxy lensing system SDP.81. The detailed morphology of the lensed multiple 
images of the z = 3.042 infrared-luminous galaxy, which is found to consist of tens of < 100- 
pc-sized star-forming clumps embedded in a ~ 2 kpc disk, are well reproduced by a lensing 
galaxy modeled by an isothermal ellipsoid with a 400 pc core. The core radius is consistent 
with that of the visible stellar light, and the mass-to-light ratio of ~ 2Af 0 Lg 1 is comparable to 
the locally measured value, suggesting that the inner 1 kpc region is dominated by luminous 
matter. The position of the predicted mass centroid is consistent to within ~ 30 mas with a 
non-thermal source detected with ALMA, which likely traces an active galactic nucleus of the 
foreground elliptical galaxy. While the black hole mass and the core radius of the elliptical 
galaxy are degenerate, a point source mass of > 3 x 10 s M 0 mimicking a supermassive black 
hole is required to explain the non-detection of a central image of the background galaxy. The 
required mass is consistent with the prediction from the well-known correlation between black 
hole mass and host velocity dispersion. Our analysis demonstrates the power of high resolution 
imaging of strong gravitational lensing for studying the innermost mass profile and the central 
supermassive black hole of distant elliptical galaxies. 

Key words: black hole physics — gravitational lensing: strong — galaxies: individual (H-ATLAS 
J090311.6+003906) — galaxies: structure — submillimeter: galaxies 
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1 Introduction 

Galaxy-galaxy lensing occurs ubiquitously in the observable 
universe. Magnified images from gravitational lensing al¬ 
low us to understand the detailed mass structure of the in¬ 
tervening galaxy as well as the background source structure 
(e.g., Koopmans 2005; Inoue & Chiba 2005; Vegetti et al. 
2012; Hezaveh 2013a; Hezaveh et al. 2013b; Hezaveh 2014). 
Furthermore, a detection or absence of a central lensed image 
offers a unique opportunity to map the innermost mass distribu¬ 
tion of the lensing galaxy and to directly measure the mass of 
the central supermassive black hole (SMBH, Mao et al. 2001; 
Winn et al. 2004; Inada et al. 2008). Precise mass modeling re¬ 
quires sensitive, extinction-free, and extremely high resolution 
imaging toward the very center of a lensing galaxy (Hezaveh 
et al. 2015). This is now possible with long baseline capa¬ 
bilities offered by the Atacama Large Millimeter/submillimeter 
Array (ALMA). Here we present a detailed modeling of a mass 
profile of the z = 0.2999 elliptical galaxy, which is the fore¬ 
ground galaxy in the galaxy-galaxy lensing system H-ATLAS 
J090311.6+003906 (SDP.81), using 30 milli-arcsec (mas) reso¬ 
lution 1-mm images obtained at ALMA. 

SDP.81 is a gravitationally lensed z = 3.042 galaxy dis¬ 
covered in the Herschel Astrophysical Terahertz Large Area 
Survey (H-ATLAS, Eales et al. 2010). The characteris¬ 
tic features (“Einstein ring”) are evident in the 2" resolu¬ 
tion interferometric maps (Negrello et al. 2010; Omont et al. 
2013; Bussmann et al. 2013) and the Hubble Space Telescope 
(HST) image (Negrello et al. 2014; Dye et al. 2014), indi¬ 
cating strong magnification. Carbon Monoxide (CO) and wa¬ 
ter (H 2 O) emission were successfully detected from the back¬ 
ground galaxy (Omont et al. 2013). The foreground elliptical 
galaxy is identified in the Sloan Digital Sky Survey (SDSS) 
as SDSS J090311.57+003906.5 (hereafter SDSS J0903) with 
a spectroscopic redshift of z = 0.2999 ± 0.0002 (Negrello et al. 
2014). We hereafter refer to the foreground and background 
galaxies as SDSS J0903 and SDP.81, respectively. Detailed dis¬ 
cussion on the spatial/kinematic structure and the star-formation 
activity in SDP.81 will be presented in a subsequent paper 
(Paper II; Hatsukade et al. 2015). 

Throughout this paper, we assume a flat universe with D m = 
0.26, Da = 0.74, and Hq — 72 km s -1 Mpc -1 . The angular 
scale of 1" corresponds to 4.36 kpc at the lens redshift z\ = 
0.2999 and 7.78 kpc at source redshift z s = 3.042. 

2 Data 

2.1 ALMA observations and images 

The ALMA data at wavelengths of 1.0, 1.3 and 2.0 mm 
were obtained during the long baseline campaign of Science 
Verification (ALMA Partnership 2015b). The resulting spatial 
resolution achieved at 1.0, 1.3 and 2.0 mm are 23 x 31 mas (PA 


= 15°), 30 x 39 mas (PA = 20°) and 54 x 60 mas (PA = 18°), 
with r.m.s. noise level of 10, 10, and 9 /xjy beam -1 , respec¬ 
tively. The CO (5^1), (8-7) and (10-9) images were also taken 
simultaneously in the 2.0, 1.3 and 1.0 mm bands, respectively. 
The CO images are uu-tapered to ~ 170 mas in order to increase 
the sensitivity to extended structure. Details of the observations 
are given in ALMA Partnership (2015a) and ALMA Partnership 
(2015b). 

A three-color image composed from the 1.0 and 1.3 mm 
data are shown in figure la. The quadruple images of SDP.81 
are clearly resolved, three of which are bridged by a 3"-long 
arc. The arcs consist of four bright knots connected by thin 
smooth filaments. The grainy clumps seen around the arcs are 
not coherent among the different frequency bands (figure la), 
and thus they are likely to be artifacts due to phase calibra¬ 
tion residuals. On the other hand, we find fainter, but sig¬ 
nificant knots embedded in the arcs that are consistent in all 
different waveband images. These are real features that are 
likely tracing the internal substructures of SDP.81. A central 
point source bracketed by the long and short arcs is not re¬ 
solved with the Band 7 beam and has a 1.0-2.0 mm power-law 
spectral index of a = —0.3 ± 0.7 (ALMA Partnership 2015a). 
The extrapolated flux at 20 cm is consistent with that of the 
cataloged radio source FIRST J090311.5+003906 (S 20 cm = 
1.21 ±0.16 mly beam -1 , Becker et al. 1995). The power-law 
spectral index from 1.0 mm to 20 cm is a = —0.64 ± 0.01 (see 
inset of figure lb), which is consistent with the spectral index 
measured for 1.0-2.0 mm. If the non-thermal emission arises 
from SDP.81, it should also dominate the emission from the 
arcs. However, the observed spectral index of the arc is posi¬ 
tive (typically a > 2), meaning that dust emission is the better 
model to explain the emission characteristics of the arc ALMA 
Partnership (2015a). Furthermore, no molecular lines of CO 
and H 2 O at z = 3.042 are found at the position. It is there¬ 
fore likely that the non-thermal emission arises from an Active 
Galactic Nucleus (AGN) of SDSS J0903. 

2.2 Hubble images 

We retrieve the archival HST image (proposal ID: 12194, 
Negrello et al. 2014) taken with the Wide Field Camera 3 
(WFC3) at 1.6 pm (F160W, H band). The point spread func¬ 
tion (PSF) is 151 mas in FWHM. We calibrate the coordinates 
of the WFC3 image using the catalog data from the SDSS Data 
Release 12 1 . We use 46 SDSS objects located in the WFC3 im¬ 
age to derive the astrometric fit using iraf. The r.m.s. errors of 
the fit are 85 mas in right ascension and 100 mas in declination. 
The nominal center of SDSS J0903 matches the position of the 
central non-thermal source in the ALMA images within 15 mas. 

We use the GALFIT software (Version 3.0.5, Peng et al. 

1 http://skyserver.sdss.org/dr12 



Publications of the Astronomical Society of Japan, (2015), Vol. 00, No. 0 


3 


Table 1. 

Lens model parameters. 
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The errors represent 68% confidence intervals estimated from a Markov-chain Monte Carlo simulation, (a) The central position of an isothermal ellipsoid in units of mas, 
relative to the position of the non-thermal source at 9 h 3 m Ilf 573, +0°39 / 6f 54, respectively, (b) Velocity dispersion of the isothermal ellipsoid (km s —1 ). (c) Ellipticity 
and its position angle (deg), (d) core radius of the isothermal ellipsoid (mas), (e) external shear and its position angle (deg), (f) third-order external perturbation and its 
position angle (deg), (g) multipole perturbation with (m, n) = (4, 2) and its position angle (deg). See Oguri (2010) for the definitions of the parameters. 



-1 0 1 2 0 0.1 0.2 0.3 0.4 

AR.A. (arcsec) AR.A. (aresec) 


Fig. 1 . (a) ALMA 3-color image of SDR81 (1.0, 1.15 and 1.3 mm for blue, green and red, respectively) overlaid with the Hubble WFC3/F160W (1.6 (im) image 
where the stellar light of SDSS J0903 is subtracted (contours). Two sets of counter-images of stellar peaks are indicated by filled and open stars, respectively. 
The synthesized beam size is indicated at the bottom-left corner. The origin of the image is taken at the position of a central compact non-thermal source, (b) 
The ALMA 1.0 mm image. Upper inset shows CO (5^1) spectra at the positions of the source A (upper) and E (lower). Bottom inset shows the spectral energy 
distribution of the central compact source, which is well fitted by a power-law function with a spectral index of -0.64 (solid line), suggesting the synchrotron 
emission, (c) The modeled brightness distribution on the image plane. The image is smoothed by a Gaussian with FWHM = 23 mas. The inner and outer 
ellipses represent radial and tangential critical curves, respectively, (d) The modeled brightness distribution on the source plane, which is 0.5" on a side. The 
star represents the source position of the stellar peaks denoted as filled stars in (a). The position of this panel is indicated as a dotted square in (c). The solid 
curves represent the caustics. The scale bar at the bottom-left corner shows a physical scale of 200 pc at z = 3.042. 
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Fig. 2. Correlation matrix of the lens parameters computed by a Marcov-chain Monte Carlo simulation. The probability distribution of two parameters are 
indicated by a greyscale map with contours of i<x, 2<r and 3rr. The diagonal panels show the probability distribution of each parameter in which the rest 
parameters are marginalized over. The asterisk in the x-y correlation plot marks the position of the AGN at the center of SDSS J0903. 


2010) to deblend the stellar light in SDP.81 from SDSS J0903. 
Following the procedure adopted by Negrello et al. (2014), two 
Sersic profiles are simultaneously fitted to the HST image. As 
a result, the radial profile of SDSS J0903 is decomposed into a 
brighter (//-band absolute AB magnitude of Mh = 17.56) ellip¬ 
tical with a Sersic index of n a = 4.42 and a fainter {Mh = 18.69) 
exponential disk with n a = 1.03. Although our total flux is 1.2x 
higher, this is consistent with the model presented by Negrello 
et al. (2014). The residual image shown in figure 1 clearly ex¬ 
hibits the stellar arcs which are significantly (> 100 mas) offset 
from the 1-mm ALMA arcs (see Paper II for the interpretation). 


3 Gravitational lens model 

3.1 Methods and results 

We use the GLAFIC software (Oguri 2010) 2 to reproduce the 
lensed image using an isothermal ellipsoid with a flat core. 
The choice of an isothermal ellipsoid is supported by the fact 
that the radial density profile of SDSS J0903 is consistent with 
p(r) ex r -2 (Dye et al. 2014; Rybak et al. 2015), which is the 
profile expected for an isothermal ellipsoid. The structure of 
SDP.81 is, however, too complex to be accounted for by a sin¬ 
gle isothermal ellipsoid alone. We examine the SDSS clus¬ 
ter catalog of Oguri (2014) and find a cluster of galaxies at 
s = 0.306 with the richness of ~ 21 and located ~ 2' away 


2 http://www.slac.stanford.edu/%7Eoguri/glafic/ 
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Fig. 3. Surface mass density versus projected radius of the best-fit mass 
model (blue shade), radial distribution of stellar light (dark yellow) and the 
point spread function (PSF) of the 1.6 pm HST image (dotted curve, arbitrary 
unit). The mass model smoothed by the HST PSF is shown in dashed curve. 
The surface density is normalized by the cosmological critical density for a 
source redshift z s = 3.042, E crlt = 4.3 x 10 9 M 0 kpc~ 2 . 

from SDSS J0903. Therefore we include the third-order and 
quadrupole perturbations of mass 3 , in addition to the standard 
external shear, to take into account the gravitational influence 
of the cluster and large-scale structure. 

We identify compact knots in the arcs that are common in 
the 1.0, 1.3, and 2.0 mm images and then identify two sets of 
counter-images of the knots in SDP.81. These are shown as 
open squares in figure lb (A and E). The counter-images of 
each set are extracted by visually identifying continuum inten¬ 
sity peaks which share the same/similar CO (5^1) and (8-7) 
spectral profiles at each of the peaks. Figure lb shows how the 
counter-images are identified; all of the CO (5-4) spectra at the 
positions of A have a peak at ~ 100 km s _1 while the peaks 
at E appear near —200 km s _1 and have broader line profiles 
than A, which likely correspond to the red and blue compo¬ 
nents reported by Omont et al. (2013) and ALMA Partnership 
(2015a), respectively. This method allows us to predict that the 
knots A are quadruple lensed images of a single source, whereas 
the knots E are double images and thus the (tangential) caustics 
should lie between A and E on the source plane. First, we ob¬ 
tain an initial mass model only with the isothermal ellipsoid and 
external shear using the positions of A and E. Then we use the 
initial mass model to identify groups of counter-images of knots 
found in the 1,0-mm continuum image. The knots within ~ 50 
mas from the predicted positions are identified as the counter¬ 
images, resulting in identifications of 46 images. We then refine 

3 We use glafic mass models cius3 for the third-order external perturba¬ 
tion and mpole with (m, n) = (4,2) for the quadrupole perturbation (equa¬ 
tion 13 of Oguri 2010). The central positions of the perturbation compo¬ 
nents are fixed to that of the isothermal ellipsoid. 


the mass model by adding the higher-order perturbations (the 
number of independent parameters is 12 ) and fitting all of the 
identified multiple images. In the model calculation, we also 
consider two sets of double images identified with the 1.6-pm 
HST image, which is indicated by star symbols in figure la. The 
number of degrees of freedom (DOF) is thus IVdof = 88. The 
positional uncertainty of the knots used for optimization is set to 
20 mas, while the uncertainties of the 10 knots which are close 
to the critical curve and/or possibly suffering from blending are 
set to 100-200 mas. The uncertainties of the four HST images 
are set to 130-200 mas. 

The best-fit parameters of the mass model 4 are listed in 
table 1, where the la confidence intervals are computed by 
Marcov-chain Monte Carlo simulations as shown in figure 2. 
Note that the core radius (r core ) and velocity dispersion (a v ) 
of the isothermal ellipsoid are highly degenerate. The central 
position of the isothermal ellipsoid is not only consistent with 
the brightness peak of the stellar distribution within the lcr as¬ 
trometric uncertainty of HST (~ 100 mas), but also with the 
location of the candidate AGN, i.e., ( x , y) = (0, 0), to within 
~ 30 mas (lcr). This suggests that the AGN is indeed located 
at the bottom (innermost < 100-pc) of the potential well of 
SDSS J0903. The best-fit parameters of the isothermal ellip¬ 
soid and external shear are in good agreement with those inde¬ 
pendently obtained using the same ALMA images (Rybak et al. 
2015; Dye et al. 2015) 5 and the HST images (Dye et al. 2014) 
at lcr level. 

The crosses in figure lb mark the positions of the counter¬ 
images predicted by this model, and we find that the model is 
qualitatively consistent with the peak positions of the 1.0-mm 
brightness distribution (\ 2 = 3.45). Figures lc and Id show 
the predicted brightness distributions on the image and source 
planes, respectively. The observed images are reproduced by 35 
< 100-pc-sized knots embedded in an extended (~ 2 kpc) disk. 
While the uncertainty is relatively large, the total magnification 
factor is expected to be p g ~ 22. This is consistent with the 
value obtained by Rybak et al. (2015) (p g = 17.6 ±0.4 for total, 
p g = 25.2 ± 2.6 for a central star-forming disk), while it is ~ 
2x higher than the value previously obtained with the SMA 
(p g ~ 11, Bussmann et al. 2013). This is likely because, as 
suggested by Rybak et al. (2015), the ALMA image allows us to 
identify more compact structures close to the caustics, resulting 
in higher magnifications than those obtained previously. 

3.2 Mass profile of SDSS J0903 

In figure 3, we show the surface mass profile of the gravitational 
lens model normalized by a cosmological critical density for a 
source redshift z s = 3.042 (i.e., convergence /t(r)) and a radial 

4 The glafic model of SDP.81 is available at httpV/www.ioa.s.u- 
tokyo.ac.jp/%7Eytamura/SDP81/. 

5 Note that the ellipticity q used in Rybak et al. (2015) is related by e = l - q. 
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SMBH mass M BH (M = ) 

Fig. 4. The flux density ratio of central counter-images relative to the all 
counter-images, as a function of mass of a supermassive black hole of 
SDSS J0903 (Solid line with circles). The source-plane components are 
the same as shown in figure Id. The horizontal dotted line shows the 30- 
upper limit of the ALMA 1.0-mm image on the flux ratio. The dashed line 
with squares represents x 2 of model optimizations. 


profile of the 1.6-pm stellar light obtained by the HSTAVFC3 
F160W imaging observations. The critical density is E cr i t = 
4.3 x 10 9 M a kpc -2 , and the stellar surface mass is scaled using 
a fiducial mass-to-light ratio of 2 M 0 Lq 1 at A re st — 1 pro 
(e.g.. Bell et al. 2003), which is a typical value for local massive 
ellipticals. Because the stellar profile is lowpass-filtered by the 
PSF of HST (151 mas, see the dotted curve), the profile is only 
reliable at r > 0.1". It is evident that there is a knee in the stellar 
profile at a characteristic radius of ~ 0.15". 

In figure 3, we also show the lensing mass profile smoothed 
by the HST PSF. The knee position and the amplitude of the 
profile are in agreement with the lens model. Although the de¬ 
rived core (t core = 48t“ 9 mas, which corresponds to a core 
diameter of ~400 pc) is slightly smaller than the stellar core 
(r ~ 150 mas), the two core radii are consistent to within the 1 <t 
uncertainty. Furthermore, the position angle of the isothermal 
ellipsoid (8 e = 11 ± 13°) is consistent with that of the stellar 
profile (10 ± 1°, Dye et al. 2014), suggesting that the bary- 
onic matter, traced by the stellar light, dominates the mass in 
the inner 1-kpc region of the galaxy. At larger radii (> 3 kpc), 
however, the stellar light falls below the model surface density 
(k(t) oc r -1 ), which can be interpreted as the presence of dark 
matter in the outer parts of SDSS 10903. 


4 Discussions and conclusions 

The non-detection of the central ‘odd’ image of SDP.81 in the 
deep ALMA images rules out a top-flat plateau in the gravita¬ 
tional potential, but instead requires the presence of a steep (or 
cuspy) radial profile in the innermost potential distribution. A 
SMBH can make the central potential more cuspy even if the 
stellar mass distribution is top-flat as suggested from the WFC3 


radial profile. As a consequence, this can erase the central im¬ 
age if the black hole is massive enough to fully perturb the core 
potential. Since there is a clear sign of the presence of an AGN 
at the center of SDSS J0903, we add a point source mimick¬ 
ing the SMBH to the lens model in order to investigate how the 
SMBH affects the flux of the central image. A point mass with 
Mbh = 10' to 10 9 M 0 is placed at the position of the AGN and 
r core is fixed to the stellar core radius of 150 mas, while the rest 
are treated as free parameters (A'dof = 89). 

Figure 4 shows S ce nt/Stot (i.e. the ratio between the flux 
from the central image to all of the counter-images) as a func¬ 
tion of Mbh- The fit to the model is generally good (x 2 — 3.7) 
in the range of 10' < Mbh/M 0 < 10 9 . The ratio is insensi¬ 
tive to the black hole mass if it is less than 1 x 10 8 Mg, but 
the central image rapidly dims as the black hole mass is higher 
than 2 x 1O 8 M 0 . The horizontal dotted line shows the 3 a up¬ 
per limit of the 1.0-mm image (30 ply) divided by the total flux 
(37.7± 1.5 mly, ALMA Partnership 2015a), suggesting that the 
ALMA non-detection places a constraint to the black hole mass 
of Mbh > 3 x 10 8 M 0 . Note that if we adopt r co re = 48 mas, 
the limit is slightly relaxed to Mbh > 1 x 10 8 M 0 . 

A black hole mass can independently be estimated from 
the well-known correlation between black hole mass and the 
velocity dispersion of the bulge (e.g., Ferrarese & Merritt 
2000; Gebhardt et al. 2000). The black hole estimated from the 
velocity dispersion of the lens (cr„ = 265j)4 5 km s -1 , Table 1) 
and using the local correlation (Kormendy & Ho 2013) is 
Mbh = (1 ± 0.1) x 10 9 M 0 , which is consistent with the lower 
limit derived above. It will difficult to detect the central im¬ 
age of SDP.81 even in future ALMA observations if Mbh ~ 
1 x 10 9 M 0 . We demonstrate in this work that the detailed grav¬ 
itational lens modeling using sensitive, high-resolution imaging 
with ALMA is capable of constraining the central < 1 kpc mass 
profile and the black hole mass, which will help us understand 
the co-evolution of the SMBHs and the host galaxies through, 
for example, SMBH scouring (Hezaveh et al. 2015). 

Finally, we note that a similar conclusion on the mass of the 
SMBH (log Mbh/M 0 > 8.4) was made by a recent analysis by 
Wong et al. (2015), where they used an independent code glee 
(Suyu & Halkola 2010) and different sets of multiple images 
from the ALMA images tapered to ~ 170 mas to construct a lens 
model. Nevertheless, we prefer to use the full resolution ALMA 
images (up to 23 mas), which allows us to (1) accurately model 
higher-order perturbations, and to (2) robustly identify multiple 
features in the images and construct a reliable lens model. 


Acknowledgments 

We thank the anonymous referee for fruitful comments. This work is sup¬ 
ported by KAKENHI (No. 25103503, 26800093 and 20647268). MO 
is supported in part by World Premier International Research Center 
Initiative (WPI Initiative), MEXT, Japan. DI is supported by the 2015 







Publications of the Astronomical Society of Japan, (2015), Vol. 00, No. 0 


7 


Inamori Research Grants Program. This paper makes use of the follow¬ 
ing ALMA data: ADS/JAO.ALMA #2011.0.00016.SV. ALMA is a part¬ 
nership of ESO (representing its member states), NSF (USA) and NINS 
(Japan), together with NRC (Canada), NSC, ASIAA (Taiwan) and KASI 
(Republic of Korea), in cooperation with the Republic of Chile. The Joint 
ALMA Observatory is operated by ESO, NAOJ and NRAO. 


References 

ALMA Partnership, Vlahakis, C., Hunter, T. R., et al. 2015, ApJL, in 
press (arXiv: 1503.02652) 

ALMA Partnership, Fomalont, E. et al. 2015, ApJL, submitted 
Becker, R. H„ White, R. L„ & Helfand, D. J. 1995, ApJ, 450, 559 
Bell, E. F.. McIntosh, D. H„ Katz, N„ & Weinberg, M. D. 2003, ApJS, 
149, 289 

Bussmann, R. S., Perez-Fournon, I., Amber, S., et al. 2013, ApJ, 779, 25 
Dye, S., Negrello, M„ Hopwood, R„ et al. 2014, MNRAS, 440, 2013 
Dye, S., Furlanetto, C., Swinbank, A. M., et al. 2015, MNRAS, submitted 
(arXiv: 1503.08720) 

Eales, S., Dunne, L., Clements, D., et al. 2010, PASP, 122, 499 
Ferrarese, L., & Merritt, D. 2000, ApJL, 539, L9 
Gebhardt, K., Bender, R., Bower, G., et al. 2000, ApJL, 539, L13 
Hatsukade, B., Tamura, Y., Iono, D., et al. 2015, PASJ, submitted 
(arXiv: 1503.07997) 

Hezaveh, Y. D„ Dalai, N„ Holder, G„ et al. 2013, ApJ, 767, 9 
Hezaveh, Y. D., Marrone, D. P., Fassnacht, C. D., et al. 2013. ApJ, 767, 
132 

Hezaveh, Y. D. 2014, ApJL, 791, L41 

Hezaveh, Y. D„ Marshall, P. J„ & Blandford, R. D. 2015, ApJL, 799, L22 
Inada, N., Oguri, M., Falco, E. E., Broadhurst, T. J., Ofek, E. O., 
Kochanek, C. S., Sharon, E., & Smith, G. P. 2008, PASJ, 60, L27 
Inoue, K. T„ & Chiba, M. 2005, ApJ, 633, 23 
Koopmans, L. V. E. 2005, MNRAS, 363, 1136 
Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511 
Mao, S., Witt, H. J., & Koopmans, L. V. E. 2001, MNRAS, 323, 301 
Negrello, M., Hopwood, R., De Zotti, G., et al. 2010, Science, 330, 800 
Negrello, M„ Hopwood, R„ Dye, S„ et al. 2014, MNRAS, 440, 1999 
Oguri, M. 2010, PASJ, 62, 1017 
Oguri, M. 2014, MNRAS, 444, 147 
Omont, A., Yang, C., Cox, P„ et al. 2013, A&A, 551, AA115 
Peng, C. Y., Ho, L. C., Impey, C. D„ & Rix, H.-W. 2010, AJ, 139, 2097 
Rybak, M., McKean, J. P., Vegetti, S., Andreani, P., & White, S. D. M. 

2015, MNRAS, submitted (arXiv: 1503.02025) 

Suyu, S. H„ & Halkola, A. 2010, A&A, 524, AA94 
Vegetti, S„ Lagattuta, D. J., McKean, J. P.. et al. 2012, Nature, 481, 341 
Winn, J. N., Rusin, D., & Kochanek, C. S. 2004, Nature, 427, 613 
Wong, K. C., Suyu, S. H., & Matsushita, S. 2015, MNRAS, submitted 
(arXiv: 1503.05558) 



